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Abstract. We report the results of XMM-Newton and BeppoSAX observations of the radio and X-ray emitting 
£^ ■ star LS I +61 303, likely associated with the gamma-ray source 2CG 135+01 and recently detected also at TeV 

energies. The data include a long XMM-Newton pointing carried out in January 2005, which provides the deepest 
, look ever obtained for this object in the 0.3-12 keV range. During this observation the source flux decreased from 

^ ■ a high level of ~13xl0 -12 erg cm -2 s _1 to 4xl0 -12 erg cm -2 s _1 within 2-3 hours. This flux range is the same 

seen in shorter and less sensitive observations carried out in the past, but the new data show for the first time 
that transitions between the two levels can occur on short time scales. The flux decrease was accompanied by a 
significant softening of the spectrum, which is well described by a power law with photon index changing from 
1.62+0.01 to 1.83+0.01. A correlation between hardness and intensity is also found when comparing different short 
observations spanning almost 10 years and covering various orbital phases. LS I +61 303 was detected in the 15-70 
keV range with the PDS instrument in one of the BeppoSAX observations, providing evidence for variability also 
in the hard X-ray range. The X-ray spectra, discussed in the context of multiwavelength observations, place some 
' interesting constraints on the properties and location of the high-energy emitting region. 
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■ 1. Introduction The hard power law X-ray spectrum extending up to 

, 10 keV without breaks (photon index of ~1.7, absorbing 

The peculiar radio source GT 0236+610, unambiguously co i umn density, N H , of -5xl0 21 cm" 2 ) is clearly incon- 



O 
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associated with the B0 Ve star LS I +61 303, is unique in sistent with low temperature plasma emission, thus ex- 
. its high variability and periodic radio emission (Gregory c i uding the Be-star as a main contribution to the X-ray 
& Taylor 1978). The radio outbursts show a periodicity 



emission. 



of about 26.5 days (Taylor & Gregory 1982; Gregory et 
al., 1999) and a further modulation of both the outburst The 26 - 5 da y s periodicity, reflecting the orbital mo- 
phase and outburst peak flux with a period of -1600 days tion of thc system (Gregory & Taylor 1978, Taylor & 
(Gregory et al. 1999), also displayed in the H Q emission Gregory 1982), has also been observed in the optical 
line (Zamanovet al 1999) band (Hutchings & Crampton 1981; Mendelson & Mazeh 

. „ . , v 1989), in the infrared (Paredes et al. 1994), in soft X- 

A taint X-ray counterpart (F 2 _iokoV, ' m , , , , . „ . . ,. 

„ 1 „_io _t _in • i ;-n i i • rays Paredes et al. 1997 and m the H„ emission line 

6x10 erg cm s 1 ) was identified by Bignami ,^ K et I 1999) 

et al. (1981) with the Einstein satellite, and later ob- ^ ' 

served with ROSAT (Goldoni & Mereghetti 1995; Taylor Spectral line observations of the radio source give a dis- 

et al. 1996), ASCA (Leahy et al. 1997), RXTE (Harrison tance of 2.0+0.2 kpc (Frail & Hjellming 1991), implying 

et al. 2000; Greiner & Rau 2001). an X-ray luminosity of 10 33 erg s -1 . No periodic pulsa- 
tions have been detected in its X-ray and radio emission. 

Send offprint requests to: L.Sidoli (sidoli@iasf-milano.inaf.it) The low X " ra y luminosity and the lack of iron line emis- 

* Based on obervations obtained with XMM-Newton, an sion indicate that LS I +61 303 is not a classical accret- 

ESA science mission with instruments and contributions di- ing X-ray pulsar. Although this is not surprising in view 

rectly funded by ESA Member States and NASA. of the strong radio emission (never seen from an X-ray 
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pulsar, see e.g. Fender 2001), the low X-ray luminosity is 
remarkable in a source that can only be explained with 
the presence of a compact object. 

LS I +61 303 is very likely associated with one of 
the most interesting unidentified gamma-ray sources near 
the plane of the Galaxy, the COS-B source 2CG 135+01 
(Hermsen etal. 1977). Although the association between 
2CG 135+01 and LS I +61 303 was initially weakened by 
the presence of another plausible counterpart in the COS- 
ES error region (the quasar QSO 0241+622), the reduced 
error box determined by EGRET (Kniffen et al. 1997) is 
compatible only with the position of LS I +61 303. 

A recent detection of variable and likely periodic very 
high energy gamma-ray emission above 100 GeV has been 
reported by the MAGIC collaboration (Albert et al. 2006). 

The strong radio outbursts suggest the presence of 
a compact star, although there is no direct evidence for 
the existence of a neutron star (no pulsations nor X-ray 
bursts) . Two main classes of models have been proposed 
for LS I +61 303, involving a highly eccentric binary sys- 
tem with an orbital period of 26.5 days, composed of a 
compact star and the Be companion. 

The first class of models suggests that the radio out- 
bursts are produced by streams of relativistic particles 
originating in episodes of super-Eddington accretion onto 
the compact object (e.g., Taylor et al., 1992). The X- 
ray luminosity of LS I +61 303 is orders of magnitude 
lower than the Eddington limit. The association with the 
gamma-ray source 2CG 135+01 suggests that the bulk of 
the energy output is shifted from X-ray to 7-ray wave- 
lengths, but the actual mechanism responsible for this is 
not well understood in the context of the supercritical ac- 
cretion model. 

The second class of models assumes that LS I +61 303 
contains a non-accreting young rapidly rotating neutron 
star (e.g., Maraschi & Treves, 1981; Tavani, 1994). The 
radio outbursts are produced by energetic electrons ac- 
celerated in the shock boundary between the relativistic 
wind of the young pulsar and the wind of the Be star. This 
scenario is supported by the shock-powered emission ob- 
served by ASCA and CGRO near the periastron in a pos- 
sibly similar system, the binary PSR B1259-63 (Tavani et 
al., 1995), which is composed by a radio pulsar (with a pe- 
riod of 47 ms) orbiting a Be star (orbital period of 3.4 yr). 
In this case, the high-energy emission allows a sensible di- 
agnostics of the shock emission region where the pulsar 
wind interacts with the circumstellar material originating 
from the surface of the massive star (Tavani, 1994). The 
modulation of the radio emission might be due to the time 
variable geometry of a 'pulsar cavity' as a function of the 
orbital phase. 

A third way of producing X-rays has been suggested by 
Campana et al. (1995), as accretion onto the pulsar mag- 
netosphere when the stellar wind ram pressure exceeds 
the pulsar wind pressure, but is not large enough to pene- 
trate the pulsar magnetospheric boundary. This situation 
can occur in LS I +61 303 for reasonable values of the 
neutron star magnetic field and spin period. 



An object similar to LS I +61 303 has been discov- 
ered, the microquasar LS 5039 (Paredes et al. 2000, 2005) 
which is subluminous in the X-ray range (even more than 
LS I +61 303) and, if the proposed association with the 
EGRET source 3EG J1824-1514 is correct, also shows the 
same puzzling behavior, having L 7 > Lx- 

Therefore, LS I +61 303 and LS 5039 could be the first 
two examples of a new class of X-ray binaries with pow- 
erful 7-ray emission (see also Aharonian et al. 2005 for a 
HESS detection of LS 5039 at energies above 250 GeV). In 
the case of LS 5039, an ejection process, probably fed by 
an accretion disk, is clearly proved by a VLBA map which 
shows bipolar jets emerging from a central core. Its high 
L 7 is tentatively explained by inverse Compton scattering. 
In LS I +61 303, VLBI observations (at 5 GHz) show a 
possible jet-like elongation at milliarcseconds scales, sug- 
gesting the presence of a one-sided collimated radio jet, 
as found in several other X-ray binaries. A lower limit of 
0.4 c for the intrinsic velocity of the radio jet has been 
derived (Massi et al. 2001). 

We report here the results of a deep XMM-Newton ob- 
servation performed in January 2005, together with spec- 
tral analysis of archival BeppoSAX observations never re- 
ported in literature, and for completeness we studied also 
other five short archival XMM-Newton observations re- 
cently analysed by Chernyakova et al. (2006). 

2. Observations and data reduction 

2.1. XMM-Newton observations 

The XMM-Newton Observatory includes three 1500 cm 2 
X-ray telescopes each with an European Photon Imaging 
Camera (EPIC) at the focus. Two of the EPIC imaging 
spectrometers use MOS CCDs (Turner et al. I2001fl and 
one uses a pn CCD (Struder et al. 12001(1 . Behind two of 
the telescopes there are Reflection Grating Spectrometers 
(RGS, 0.35-2 keV; den Herder et al. 2001). LS I +61 303 
was observed on 2005 January 27-28 for a net exposure 
of 48.7 ks. Data were processed using version 6.5 of the 
XMM-Newton Science Analysis Software (SAS). Known 
hot, or flickering, pixels and electronic noise were rejected 
using the SAS. A further cleaning was necessary because 
of the presence of a soft proton flare at the beginning of 
the observation, reducing the net good exposure times to 
42.6 ks and 30.3 ks, respectively for the EPIC MOS and 
the pn. EPIC MOS used Imaging Prime Partial Mode, 
while EPIC pn the Small Window Mode. All EPIC cam- 
eras used the medium thickness filter. 

Cleaned MOS (with pattern selection from to 12) and 
pn (patterns from to 4) events files have been extracted, 
and used for the subsequent analysis with XMMSELECT. 
The source spectra have been extracted from a circular 
region with 40" radius, for the MOS1, MOS2 and pn sep- 
arately. Background spectra have been always extracted 
from source free regions of the same observation. When 
time selected spectra have been considered, the same time 
selection has been applied also to the background. Spectra 
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have always been rebinned such that at least 20 counts per 
bin were present and such that the energy resolution was 
not over-sampled by more than a factor 3. The response 
matrices have been generated using the rmfgen and AR- 
fgen tools implemented in the SAS 6.5. Free relative nor- 
malizations between MOS1, MOS2 and pn instruments 
were included in the spectral fitting. All spectral uncer- 
tainties and upper-limits are given at 90% confidence for 
one interesting parameter. 

We also analysed other 5 much shorter XMM-Newton 
observations, performed in 2002. Table^gives a summary 
of the observations and the corresponding orbital {9\) and 
super-orbital phases (O2 )■ The mean error on the orbital 
phase is ±0.03. Here we adopt as time of the zero phase 
JD =2,443,366.775, together with the updated estimates 
of the orbital (Pi=26.4960±0.0028 days) and of the super- 
orbital period (a periodic modulation of the phase and 
amplitude of the radio outbursts, with a period, P2, of 
1667±8 days) reported in Gregory (2002). 

2.2. BeppoSAX observations 

LS I +61 303 was observed twice with BeppoSAX, on 1997 
September 22 at 19:43 and on 1997 September 26 at 14:55 
(UTC), with an on-source time of 12 ks and 8.6 ks, re- 
spectively We report the results obtained with the Low- 
Energy Concentrator Spectrometer (LECS; 0.1-10 keV; 
Parmar et al. I1997fl . the Medium-Energy Concentrator 
Spectrometer (MECS; 1.8-10 keV; Boella et al. ITMTjl 
and the Phoswich Detection System (PDS; 15-300 keV; 
Frontera et al. H997f) . 

Counts for the spectral analysis have been extracted 
from circular regions with standard radii (4' for the MECS 
and 8' for the LECS). The LECS and MECS spectra were 
rebinned to oversample the full width half maximum of 
the energy resolution by a factor 3 and to have a minimum 
of 20 counts per bin. Background spectra were extracted 
from annular regions around the source, in the same obser- 
vations. Response matrices appropriate for the sizes of the 
extraction regions were used. In the spectral analysis we 
fitted together LECS (selected in the best calibrated en- 
ergy range 0.1-4.0 keV) and MECS spectra (1.8-10 keV), 
with free relative normalizations between the instruments. 

The PDS consists of four independent non-imaging 
units arranged in pairs each having a separate collimator. 
Each collimator was alternatively rocked on- and 210' off- 
source every 96 s during the observations. After a proper 
cleaning of the events files aimed at excluding the "spikes" 
due to single particles hits, we extracted net spectra from 
the two observations. Inspection and comparison of the 
background spectra show that there are no known con- 
taminating sources in the background region. Since the 
only nearest possible high energy contaminating source 
is QSO 0241+622, which is located just outside the PDS 
field of view, at 1.4°, we are confident that the high energy 
emission is produced by LS I +61 303. 
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3. Results 

3.1. The XMM-Newton observations 

The EPIC pn background subtracted lightcurves in two 
energy ranges (0.3-2 keV and 2-12 keV) are reported 
in Fig. 2] together with their hardness ratio (HR). Both 
lightcurves show a smooth modulation, ranging from 
a maximum rate around 1 counts s _1 down to about 
0.2 counts s , in few hours. The hardness ratio shows 
a rather constant value of 1 for the first 15,000 s, and 
then a sharp decrease, where it falls down to 0.8 in less 
than 1000 seconds. After this decrease, the hardness ra- 
tio is smoothly modulated, displaying a wave-like shape, 
with two maxima (around 30,000 s and 46,000 s in Fig.Q 
and one minimum. A plot of HR versus the source inten- 
sity (Fig. |5J indicates that the source is harder when it is 
brighter. 



Since HR is variable, we extracted two spectra, one 
corresponding to the first 14,000 s, when HR was con- 
stant and the second one for the remaining part of the 
observation. In both cases, the continua appear feature- 
less, and an absorbed power-law model gives the best-fit. 
There is no evidence for a break up to 12 keV, and a 
bremsstrahlung model is unacceptable, with a reduced x 2 
larger than 1.4. An additional component (a blackbody) 
added to the power-law continuum is not formally required 
by the residuals, and can contribute at most 1% of the flux 
observed in the energy range 0.5-10 keV. Adopting a sin- 
gle absorbed power-law, there is no evidence for column 
density variations between the two time-selected spectra, 
thus we fixed it at 5xl0 21 cm -2 . The best-fit spectral re- 
sults are reported in Table [3 and shown in Fig 02 

The RGS spectra extracted from the whole observation 
appear featureless, with a count rate of 0.019 counts s _1 
and 0.025 counts s _1 in the RGS1 and RGS2, respec- 
tively (0.3-2 keV). The combined RGS1 plus RGS2 
spectra are well fitted (x 2 /dof=110. 7/127) with an ab- 
sorbed powerlaw with an absorbing column density of 
(6±l)xl0 21 cm -2 , and a poorly constrained photon in- 
dex of 1.9+0.6. The statistics is too low to allow a 
time/hardness selected RGS spectral analysis, so we will 
not discuss them further. 

A search for pulsations resulted in negative results. 

We have reduced and analysed in a similar way also 
the five shorter XMM-Newton observations. Within each 
observation there is no evidence for significant changes in 
hardness ratio, so we extracted spectra from the whole 
exposures (after cleaning for possible contamination from 
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Table 1. Summary of the XMM-Newton observations reported here. "Exp" is the nominal exposure time of a single 
MOS camera. 0\ and 02 are the orbital and super-orbital phases (Gregory 2002). 



OBS ID. 


Start time 


End time 


Exp. 


6>i 


02 




(dy mon yr hr:mn) 


(dy mon yr hr:mn) 


(ks) 






00 


27 Jan 2005 17:41 


28 Jan 2005 07:37 


48.7 


0.61 


0.02 


01 


10 Feb 2002 11:26 


10 Feb 2002 12:50 


6.4 


0.76 


0.37 


02 


05 Feb 2002 06:24 


05 Feb 2002 03:30 


6.4 


0.55 


0.37 


03 


21 Feb 2002 16:02 


21 Feb 2002 17:26 


7.5 


0.18 


0.38 


04 


17 Feb 2002 04:40 


17 Feb 2002 06:04 


6.4 


0.01 


0.37 


05 


16 Sep 2002 02:29 


16 Sep 2002 04:09 


6.4 


0.97 


0.50 




Time (s) 

Fig. 1. Background subtracted EPIC pn lightcurves in two energy ranges (soft=0.3-2 keV; hard=2-12 keV) together 
with their hardness ratio (HR=2-12 keV / 0.3-2 keV). Time is in seconds starting from 2005, January 27, at 17:56:26. 
Bin time is 1024 s. 



proton flares). The spectral results obtained fitting the X- 
ray emission with an absorbed power-law (and absorbing 
column density fixed at 5xl0 21 cm -2 ), are reported in 
Table El 

3.2. BeppoSAX observations 

The background subtracted lightcurves of LS I +61 303 
during the two BeppoSAX observations are shown in 
Fig. The source was on average ~2.5 times brighter 
during the second observation. Within each BeppoSAX 
observation the source displays some evidence for flux vari- 
ability, especially at the beginning of the first observation, 
where a sort of "flare" is present, with the count rate in- 
creasing from ~0.05 counts s -1 to ^0.23 counts s _1 during 
the first 1500 s, without evidence, within the uncertainties, 
for a change in hardness ratio. 



An absorbed power-law is the best-fit to both spec- 
tra (LECS+MECS2+MECS3), with the parameters re- 
ported in Tab. We fixed the absorbing column density 
to 5xl0 21 cm -2 . A bremsstrahlung model is not a good 
fit to the spectra, resulting in a temperature higher than 
18 keV, well above the MECS spectral range. A single 
absorbed blackbody is unacceptable, with reduced x 2 >2. 

The analysis of the PDS data resulted in a detection 
at high energy (15-70 keV) only in the second observa- 
tion, which indeed catched the source at a higher flux 
level, almost three times higher than in the first one. 
The PDS count rates are 0.098+0.048 counts s" 1 and 
0.31+0.07 counts s , respectively for the first and sec- 
ond observation (15-70 keV). Thus, only the second can 
be considered a detection (at 4.4 a). This PDS rate trans- 
lates into a 15-70 keV flux of about 3x 10 -11 erg cm" 2 s" 1 , 
adopting the best-fit reported in Table [3 
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Fig. 2. Hardness ratio (HR) versus Intensity (count rate in the total energy range 0.3-12 keV) plot for the longest 
XMM-Newton observation (EPIC pn), with a bin time of 1024 s. The hardness ratio is defined as ratio between net 
counts in the hard range 2-12 keV and net counts in the soft energy range 0.3-2 keV. 



Table 2. Best-fit results of the overall spectral anal- 
ysis. XOOfirst and XOOsecond indicate the two time- 
selected spectra (MOSl+MOS2+pn, 0.3-12 keV) of the 
longest XMM-Newton observation, X01 — X05 mark the 
five short XMM-Newton observations, while 501 and 
S*02 indicate the results from the two BeppoSAX point- 
ings (LECS+MECS2+MECS3; 0.1-10 keV). The best-fit 
model is and absorbed power-law, with the absorbing col- 
umn density fixed at 5xl0 21 cm~ 2 . Fluxes are corrected 
for the absorption and are in units of 10 -12 erg cm~ 2 s _1 
in the energy range 2-10 keV. 



Spectrum 


Photon 
Index 


Flux 


X 2 /dof 


XOOfirst 


1.62±0.01 


12.9±0.1 


680.7/596 


XOOsecond 


1.83±0.01 


4.0±0.1 


617.5/533 


X01 


1.53±0.02 


13.5±0.1 


494.6/462 


X02 


1.60±0.02 


13.9±0.1 


489.1/442 


X03 


1.57±0.07 


5.0±0.1 


140.1/174 


X04 


1.74±0.05 


6.8±0.1 


182.2/178 


X05 


1.52±0.02 


13.6±0.1 


497.6/453 


S01 


i-68±8:ig 


4.8±0.1 


20.0/31 


S02 


1.56±0.09 


14.0±0.1 


60.8/62 



4. Discussion and Conclusions 

We have reported here the deepest X-ray observation of 
LS I +61 303 ever carried out. A systematic analysis of 
other five XMM-Newton and two BeppoSAX observations 



has also been performed, in order to investigate possible 
changes of the spectral parameters along the orbital phase 
(see Figs. and EJl and to get a global picture at X-rays. 

During the longest XMM-Newton observation for the 
first time we found evidence for a rapid (on timescales of 
few hours) change in flux and hardness ratio: the source 
flux decreased by a factor of ~3, together with a drop in 
hardness ratio (see Fig. ^| within about 1000 s. 

This kind of variability (the source is harder when is 
brighter) is also visible when comparing different obser- 
vations performed in different times and with different 
satellites (FigEJ). Another interesting variability feature 
is the presence of a sort of "flaring" behaviour at the be- 
ginning of one of the BeppoSAX observations, but with 
no evidence for simultaneous hardness variations, mainly 
because of the low statistics. 

The evolution of the source X-ray flux with the orbital 
phase shows X-ray emission along all the phases, although 
with a difference of a factor 3 in the flux level between 
two states, with a "high" state preferentially found in the 
phase range 0.4-1.0. We found no obvious correlation of 
any of the X-ray spectral parameters with the superorbital 
period. 

The analysis of archival BeppoSAX observations al- 
lowed the first detection of the source with the PDS in- 
strument (in the energy range 15-70 keV), with evidence 
for a change in flux between a high and low state also in 
the hard energy band. 
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Fig. 3. Left panels: Two LS I +61 303 time selected counts spectra (MOSl+MOS2+pn), together with the residuals 
in units of standard deviations when fitted with an absorbed powerlaw. From top to bottom, the first and the second 
spectra are shown. Right panels: The two corresponding time selected photon spectra. 



The observations reported here demonstrate that the 
source is variable both at soft and hard X-rays, also at 
short timescales. Thus, a detailed study of the overall spec- 
trum should be in principle performed with simultaneous 
observations at all wavelengths. We try here to discuss 
the broad band energy spectrum distinguishing at least 
two source states, a high and a low state, in order to get 
informations on the source energetics. 

Figured shows °ur XMM-Newton and BeppoSAX re- 
sults in the multiwavelength spectral energy distribution 
(SED) of LS I +61 303. For radio, XMM-Newton and 
BeppoSAX data we show both high and low state spec- 
tra. CGRO and MAGIC observations show low signifi- 
cance flux variations (peaking in the phase interval 0.4- 
0.7) correlated with those at X-ray energies. For these 
instruments we plot only average measurements of the 
datasets with positive detections. During the high state, 
the kcV-MeV spectrum can be fit with a hardening power 
law (r=1.6-1.5) up to the break in the COMPTEL energy 
range (~1-10 MeV). Then the spectrum shows a shallow 



roll-off (EGRET: L=2.1; MAGIC: r=2.6) extending up to 
TeV energies. 

Two main models to account for the SED of 
LS I +61 303 have been proposed. The "micro-quasar" 
model involves streams of relativistic particles originat- 
ing in episodes of super-Eddington accretion onto a com- 
pact star embedded in the mass outflow of the B-star 
(Taylor & Gregory 1984; Massi et al. 2004a, Massi 2004b). 
Alternatively LS I +61 303 might contain a non-accreting 
young pulsar in orbit around the mass-losing B-star 
(Maraschi & Treves 1981; Dubus 2006). In both cases, 
the high energy emission offers a diagnostics of the shock 
emission parameters in the jet or in the region where the 
hypothetical pulsar wind would interact with the circum- 
stellar material. 

The spectral index and break of the keV-MeV emis- 
sion is consistent with that of diffusive shock acceleration 
models in the "slow cooling" regime (i.e. far from energetic 
equilibrium between particle injection and cooling, see e.g. 
Chevalier 2000). The shallow roll-off in the MeV-TeV en- 
ergy range requires further model components in addition 
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Fig. 4. The BcppoSAX MECS2+MECS3 background subtracted lightcurves during the two observations (on the left 
and on the right, respectively). The bin size is 512 s. A large flux variability is visible at the beginning of the first 
BeppoSAX observation, with no evidence for simultaneous hardness ratio variations. 



to the synchtrotron emission from shock-accelerated par- 
ticles. However, the large error bars of the observations 
in this energy range can hardly constrain the several pos- 
sibilc contributions based both on leptonic mechanisms 
as Compton-scattered stellar radiation, synchrotron self- 
Compton, pair cascades (Dubus 2006; Dermer & Bottcher 
2006; Bosch- Ramon et al. 2006; Bednarek 2006), and on 
hadronic mechanisms as proton interactions, producing 
gamma-rays via neutral pion decay (Romero et al. 2003, 
2005). 

A clear signature of the SED is instead the spectral 
break in the 1-10 MeV range, possibly related to the 
maximum energy achievable from electrons in the shock 
acceleration accounting for synchrotron and Compton- 
scattered stellar radiation energy losses, the two major 
electron cooling mechanisms in the dense LS I +61 303 en- 
vironment. The maximum Lorentz factor of shocked par- 
ticle is obtained by solving the equation: 



1 



1 



^synch t, 



l 



(i) 



where t acc = £i gyr = m c c"f^/eB is the electron accelera- 
tion time to reach a Lorentz factor 7 in a magnetic field 
with mean intensity B, and £ ~ 1 because the accelera- 
tion time in Fermi processes cannot be smaller than the 
gyration time, i.e. an electron gains at most a fraction of 
its energy when executing a single gyration in first-order 
shock or second-order stochastic Fermi processes (see e.g. 
Dermer & Bottcher 2006). i syn ch = 67rm c/<7T-B 2 7 is the 
synchrotron cooling time for an electron with Lorentz fac- 
tor 7, and i comp accounts for inverse Compton energy 
losses. The solution of equation (1) is: 



7n 



67re 



CT T CB(l+t synch /t comp ) 



1.2x10 s 



(2) 



compj 



If isynch <C i C omp then the maximum synchrotron 
photon energy does not depend on the magnetic field 



B: /ifmax = heBS^^/nieC ~ 167<5/£ MeV, where S 
is the Doppler factor to be included for jet scenarios 
(1 < S < 2, Chernyakova et al., 2006). Thus, neglect- 
ing inverse Compton losses, the spectral break frequency 
would be more than one order of magnitude higher than 
that observed. 

Since photons from the Be star have a mean energy of 
2.7k B T ~ 6.5 eV (T s = 2.8 x 10 4 K), the turning point 
between Thompson and Klcin-Nishina inverse Compton 
regime is for 7kn = m c c 2 jhv ~ 8 x 10 4 . Even in case 
of severe inverse Compton energy losses, eq. 2 shows that 
7max cannot be lower than 7kn for typical magnetic fields 
values expected in microquasar jets or pulsar-stellar winds 
shocks (up to several Gauss). Thus f comp = £kn and con- 
sidering electrons with Lorentz factor 7 interacting with 
a blackbody surface radiation field (Blumcnthal & Gould 
1970): 



*KN = 7 



(A) 



2 »iSfA»(lS?) 2 [l»(( | -552ATS 



(3) 



10- 3 V[cm] 2 



where R s — 13.4i? Q is the Be star radius, r is the dis- 
tance between the star and the shock region and A - 1 
accounts for the Lorentz boosting of the stellar radiation 
field into the co-moving frame of the emission region (in 
the microquasar scenario). 

Substituting cq. 3 in eq.2 we obtain a general expres- 
sion for the maximum Lorentz factor of the electrons ac- 
counting both for synchrotron and inverse Compton losses 
in the KN regime: 



7m ax — 



'1.44 x 10 16 B[G] - 8 x 10 38 £r[cm]- 



which corresponds to a break energy of 



167 
T 



10 



25 



B[G]r[cm] : 



MeV 



(4) 



(5) 
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The observed break energy at 1-10 MeV constrains 
the value of Br 2 « 10 23 G cm 2 . Shock models associated 
to pulsar wind interacting with the circumstcllar environ- 
ments predict a standoff distance of the order of r m 10 11 
cm and magnetic fields of a few Gauss (Dubus 2006) in 
good agreement with the above constraint. 

On the other hand, in accreting models, a shock orig- 
inated in an extended jet (r w 10 14 — 10 15 cm; Massi et 
al. 2001, 2004a) would imply unrealistic low values for the 
magnetic fields, unless other cooling processes apart from 
inverse Compton on stellar radiation could play a major 
role (e.g. SSC, Gupta & Bottcher 2006). Shocks originat- 
ing in the inner jet at a distance comparable with the 
compact object-star separation (r w 10 12 cm) would im- 
ply more realistic values for the magnetic fields. 

Future deeper and simultaneous observations in the 
MeV-GeV energy range with AGILE and GLAST satel- 
lites and in VHE gamma-rays with MAGIC and HESS 
telescopes will better assess the IC parameters (and/or 
hadronic mechanisms) and thus the actual emission re- 
gion properties. 
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Fig. 5. Summary of the dependence with the orbital phase of the X-ray spectral parameters of all the observations 
analysed here (six XMM-Newton and two BeppoSAX observations) together with the ASCA results (published by 
Leahy et al. 1997) for completeness. The meaning of the symbols is as follows: thick circles represent the longest XMM- 
Newton observation, splitted into two (the "first" and the "second" spectra, with significantly different hardness ratio 
and flux); the other 5 XMM-Newton observations are marked with thick stars, the 2 BeppoSAX observations are 
indicated with thick triangles and the 2 ASCA observations (from Leahy et al. 1997) with thin crosses. Fluxes are 
corrected for the absorption and in units of 10~ 12 erg cm' 2 s _1 . 
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Fig. 6. Unabsorbed flux (2-10 keV; in units of 10~ 12 erg cm" 2 s" 1 ) versus power-law photon index, collecting all 
X-ray observations analysed here, together with the ASCA results (published by Leahy et al. 1997). The meaning of 
the symbols is the same as in Fig. [5] The source appears to be softer when it is fainter. 
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Fig. 7. Broad band LS I +61 303 spectrum from radio wavelengths to TeV energies. Radio data for high and low states 
(open triangles), IR data (open stars) and OSSE data (filled triangles) are from Strickman et al. 1998. XMM-EPIC 
data (crosses) and SAX-PDS data (filled circles) both high and low states are from this work. COMPTEL data (filled 
squares) are from van Dijk et al. 1996. EGRET data (upside-down filled triangles) are from Kniffen et al. 1997. MAGIC 
data (filled stars) are from Albert et al. 2006. 



